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ABSTRACT 

Boron abundances have been derived for seven main sequence B-type stars from 
HST STIS spectra around the B III 2066A line. In two stars, boron appears to be 
undepleted with respect to the presumed initial abundance. In one star, boron is 
detectable but it is clearly depleted. In the other four stars, boron is undetectable 
implying depletions of 1 to 2 dex. Three of these four stars are nitrogen enriched, 
but the fourth shows no enrichment of nitrogen. Only rotationally induced mixing 
predicts that boron depletions are unaccompanied by nitrogen enrichments. The 
inferred rate of boron depletion from our observations is in good agreement with 
these predictions. Other boron-depleted nitrogen-normal stars are identified from the 
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literature. Also, several boron-depleted nitrogen-rich stars are identified, and while all 
fall on the boron-nitrogen trend predicted by rotationally-induced mixing, a majority 
have nitrogen enrichments that are not uniquely explained by rotation. 

The spectra have also been used to determine iron-group (Cr, Mn, Fe, and Ni) 
abundances. The seven B-type stars have near solar iron-group abundances, as 
expected for young stars in the solar neighborhood. We have also analysed the halo 
B-type star, PG 0832-^676. We find [Fe/H] = -0.88 ±0.10, and the absence of the 
B III line gives the upper limit [B/H] < —2.5. These and other published abundances 
are used to infer the star's evolutionary status as a post-AGB star. 

Subject headings: stars: abundances, evolution, rotation 

1. Introduction 

The light trace elements lithium, beryllium, and boron (LiBeB) are at the center of 
astrophysical puzzles involving sites as diverse as the primordial fireball, interstellar or even 
intergalactic space, and stellar surfaces and interiors. This paper is concerned primarily with 
boron's role in testing models of stellar interiors and evolution. This role arises because boron 
nuclei are destroyed by warm protons, and thus even quite shallow mixing of the atmosphere with 
the interior reduces the surface abundance by bringing boron depleted material to the surface. 
Lithium and beryllium are similarly affected. Thus, a determination that the surface abundance 
of a light element is less than the star's initial abundance is an observational constraint for testing 
models of stellar interiors. 

Lithium has served to test models of cool stars at essentially all phases of evolution, from the 
pre-main sequence through to the post-AGB stage, with the atomic resonance doublet conveniently 
placed at 6707 A. Beryllium has seen similar but limited use through the Be II 3130 A resonance 
doublet. Boron has seen little use in testing mixing and other processes (e.g., mass loss) owing to 
the location of its resonance lines of B I, B II, and B III in the ultraviolet. This is unfortunate 
because boron alone is observable in hot stars. 

In the case of hot stars, boron abundances have been determined using the B II 1362 A 
resonance line by Bocsgaard & Heacox (1978) using Copernicus data, Venn, Lambert, &: Lemke 
(1996) using lUE data, and Cunha et al. (1997) from the GHRS on the Hubble Space Telescope. 
However, in stars of mid-B and earlier spectral types, boron is predominantly present as B III, 
not B II. Proffitt et al. (1999) was the first to use the B III 2066 A line to determine boron 
abundances and ^^B/^°B isotopic ratios in three B-type stars from very high S/N HST GHRS 
spectra. Subsequently, Proffitt and Qiiiglcy (2001) mined the lUE archive for high-resolution 
spectra of B-type stars, and reported boron abundances or upper limits for 44 early B-type stars. 

A principal goal of most of these studies of hot stars was to establish the present-day boron 
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abundance in order to improve our understanding of the galactic chemical evolution of boron. In 
this paper, our primary goal is different; we are searching for boron-dcplctcd stars to understand 
how their depletion arose. Only Proffitt &: Quigley (2001) have shared this goal. Herein, we shall 
describe a test of rotationally-induced mixing in rapidly rotating massive stars. Shallow mixing 
results in the destruction of boron (lithium and beryllium, also) by proton capture. Deeper mixing 
also penetrates regions where H-burning via the CN-cycle has converted carbon into nitrogen. 
Prom any deep mixing, one may thus anticipate a drop in the surface abundance of boron, followed 
by a continued decline in boron accompanied by a nitrogen enrichment and carbon depletion as 
the star ages on the main sequence. 

New stellar evolution models that include the effects of rotationally-induced mixing (Heger 
&; Langer 2000, Maeder & Meynet 2000) predict that boron is depleted during main sequence 
evolution. Fliegner, Langer, &; Venn (1996) first predicted a boron-to-nitrogen relationship 
in rotating stars and suggested it as an observational test of rotational mixing in hot stars. 
This test is unique since boron and nitrogen abundance variations, and the boron-to-nitrogen 
relationship predicted, cannot be explained by simple initial abundance variations, nor binary 
evolution with mass transfer (e.g., Wellstein et al. 2001, Wellstein 2000, discussed further below). 
Extraordinary mass loss rates for B-type main sequence stars are demanded if surface boron 
depletions are to occur through exposure of boron depleted layers. The lifetime of ^^B and ^''B 
against proton capture equals the main sequence lifetime of about 10^ years at an internal layer 
with a temperature of about 7 xlO^ K. Above this layer resides about ^1 Mq of material (see 
Fig. 1 in Fliegner, Langer, & Venn 1996, and Fig.s 2, 3, & 4 in Hcgcr & Langer 2000), and 
therefore a mass loss rate of '-^10"''' Mg/yr is required if surface boron depletions are to occur. 
Such a rate is typical of B-type supcrgiants, but at least an order of magnitude larger than derived 
for main sequence OB stars (Cassinelli et al. 1994). The new rotating stellar evolution models 
also address other long-standing problems (beyond abundance anomalies), such as the origins of 
B[e] and WNL/Ofpe (slash) stars, the distribution of blue-to-red supergiants in the HR diagram, 
and the unseen post-main sequence gap predicted in all standard stellar evolution scenarios (see 
Langer & Heger 1999). 

In this paper, we present boron abundances from the B III 2066 A line from new HST STIS 
spectra of seven Galactic main sequence B-type stars. These stars were selected from a variety of 
OB associations, and have been well-studied in the optical so that atmospheric parameters and 
surface nitrogen and carbon abundances are available in the literature. We have also used the 
spectra to determine their iron-group abundances, which are difficult to determine accurately from 
optical studies of B-type stars, due to a lack of Fe III lines and uncertainties in NLTE effects on 
available lines. 

Our observing program included an eighth star: PG 0832-1-676, a B-type star in the halo, 
whose evolutionary status has been debated for a number of years. Brown et al. (1989) deduced 

that the star was a normal Population I Bl V object at a galactoccntric distance of Re ~ 37kpc, 
and z ~ 18kpc from the galactic plane. Hambly et al. (1996) reanalysed this star and found 
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that a-elements (Mg, Al, Si, S) were underabundant by about 0.4 dex relative to HR 1886, a 
local B-type star. These authors argued that the star was not a H-burning star, but rather was 
either a blue horizontal branch or a post-AGB star. Our spectra provide not only an estimate 
of the boron abundance but also new iron-group abundances (Hambly et al. detected one Fe III 
line, which suggested [Fe/H]~ —0.5 from a differential analysis). We will use these abundances 
to further establish the star's evolutionary status. Also, as Hambly et al. noted, PG 0832+676 is 
extremely sharp-lined, and as such its ultraviolet spectrum could provide a superb template for 
line identifications and wavelength calibrations. 



2. Target Selection and Observations 

Eight B-type stars were selected for HST STIS spectroscopy near the B III 2066A hne, 
see Table |l|. Five of these targets are bright stars in OB associations. Our goal was to select 
well-studied stars with a range of carbon and nitrogen abundances, and located in OB associations 
of different ages, in order to search for trends in the boron abundances with these parameters. 
A sixth target is the bright star, HD 34078, a runaway OB star from Ori OBI. These six bright 
objects are also known variable stars (typically, P Cep variables - this seems to have no influence 
on the surface abundances, discussed below), and two are in known binary systems, see Table |2[ 
Two fainter targets are also included; BD+56° 576 in Per OBI, and PG 0832+676 in the Galactic 
halo. All of our targets have well determined atmospheric parameters, see Table ^. 

Careful planning was necessary to avoid the MAMA-NUV brightness limits in observing 
the six bright targets. Even the most narrow slit, with the high resolution E230H grating 
(R = 114,000), was rarely permitted for our observations. This combination worked only for 
HD 34078 at V = 6.0. For the other brighter objects, it was necessary to use either the lower 
resolution E230M grating (R = 30,000), or E230H with a neutral density filter, 0.2x0.05ND 
(ND=2), to obtain our spectra in the minimum number of orbits; see Table [l|. The E230M grating 
is sufficient to resolve the B III line in all of our targets; the instrumental broadening for this 
grating is 9 km s~^ (±1 km s~^ depending on the slit choice), whereas the targets have measured 
rotational velocities ranging from 1< vsini <39 km s-^ For the faint star, PG 0832+676, the 
E230M grating and 0.2x0.2 slit combination produced a good quality spectrum for the abundance 
analysis, but we were not able to take full advantage of the sharp-lined nature of PG 0832+676 as 
a UV spectral template. 

Multiple exposures were taken at slightly varying central wavelengths to maximize the S/N 
at the B III 2066A line for most objects. The exposure times and central wavelengths are listed 
in Table ||. Spectra were reduced using the STIS pipeline. Each echelle observation included 
between 34 and 52 orders. We selected the order(s) that included the B III line in each exposure, 
and combined these for maximum S/N near 2066 A. The spectra were smoothed (3-pixel boxcar 
smoothing), rectified using a low-order Legendre polynomial, and offset from vacuum wavelengths 
(observed) to air wavelengths (line list, discussed below). Sample spectra are shown in Figures |l], 
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^, and|^ The resultant peak S/N values near 2066 A are listed in Table ||. The final spectra used 
for iron-group and boron abundance determinations in this paper range from 2045 to 2078 A, with 
a dispersion ranging from ~4 (E230H) to ~16 km s~^(E230M). 

Four of the stars in this study are known (3 Cep variable stars (see Table |2|). In fact, 
HD 205021 is [3 Cep. These stars exhibit variability in their radial velocities as well as their 
luminosities. For instance, HD 205021 displays a peak-to-trough amplitude of ~40 km s~^over a 
period of 0.19 d (Hadrava & Harmanec 1996). This period and amplitude is typical of the (5 Ceph 
variables in our target list, and is long enough to display radial velocity variations over the full 
range of our observations. However, the observations were divided into very short exposures to 
avoid the MAMA brightness limits. Thus, there is little significant broadening due to the radial 
velocity variations per exposure (typically about 4 km s~^, although the longest exposure for 
HD 16582 may be broadened by as much as 9 km s~^). As an additional test, we examined the 
FWHM of our individual exposures and found insignificant differences in most cases. In the case of 
HD 216916, we were surprised to find no significant differences in the FWHM when comparing the 
262 s, 880 s, and the 1200 s exposures (although the amount of broadening due to radial velocity 
variations will also depend on the phase when the observation is taken). Finally, we shifted each 
spectrum into the stellar rest frame before coaddition to minimize smearing effects. 

3. The Abundance Analyses 

Elemental abundances have been determined from LTE spectral syntheses and ATLAS9 
model atmospheres (Kurucz 1979, 1988). Solar metallicity models were used throughout, with the 
exception of PG 0832-1-676 for which [Fe/H] = —1.0 models were adopted. LTE spectral syntheses 
were made using the program LINFOR0. 

The stellar Tefj , gravity, and projected rotational velocities (vsini) values were adopted from 
the literature, see Table ^. We note here that there is some uncertainty in the temperature scale 
for B-type main sequence stars. This is discussed further in Section 4. Other parameters were 
determined from the syntheses, i.e., macroturbulence, microturbulence (^), and radial velocity, 
see Table ^. Macroturbulence {^mo) was initially set to the instrumental broadening values, then 
increased to best fit the spectral lines. In most cases, the macroturbulence is only 3 to 5 km s~^ 
larger than the smoothed instrumental profiles, however HD 205021 and HD 16582 required 
macroturbulence values of 8 km s^^ larger. Observations for these latter two stars were made 
with E230H, which has a very sharp instrumental profile, thus much of the broadening in these 
two stars is likely due to their j3 Ceph radial velocity variations (discussed above). Another star, 
HD 34078, observed with E230H has a much lower macroturbulence and is not a known (3 Ceph 



^LINFOR was originally developed by H. Holweger, W. Steffen, and W. Steenbock at Kiel University. It has been 
upgraded and maintained by M. Lemke, with additional modifications by N. Przybilla. 
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variable. A gaussian line-spread function was adopted throughout our spectrum S3mtheses; while 
acceptable, we recognize that line spread functions for the echelle gratings are available from 
STScI, but we neglected this specific instrumental broadening component in our syntheses. 

^ was determined for each star using the Fe III lines, requiring that strong and weak lines 
yield the same iron abundance. Our ^ values are significantly lower than those in the literature 
from analyses of optical spectral lines (i.e., O II). This is not uncommon, as previously found by 
Venn et al. (1996) and Cunha et al. (1997), where ^ from UV analyses is significantly smaller than 
from optical analyses of the same star. This probably reflects differences in the structure of the 
upper atmosphere where the UV spectrum forms. Regardless, the value of ^ has little impact on 
the derived B III 2065.8 A line abundance (see below). 

Only the BD-|-56° 576 vsini value was determined from syntheses since the upper-limit in 
the literature was clearly too large. BD-|-56° 576 and PG 0832+676 were observed with the same 
instrumental set-up, thus the broadening parameters from the PG 0832+676 analysis were adopted 
for the BD+56° 576 analysis. PG 0832+676 is a known sharp-lined object (fsim=l km s~^), thus 
any remaining broadening required for BD+56° 576 was attributed to its projected rotation rate. 

3.1. Line List 

Our spectral line list and atomic data originated from the Kurucz (1988; CD-18) line list, 
including all lines in the the iron-group, light elements, and heavy elements lists, up to barium, and 
through the flfth ionization state from titanium onwards. We updated this line list by including 
the new wavelengths for eight Fe III lines reported by Proffitt et al. (1999) from FTS laboratory 
measurements. We also updated the atomic data from Kurucz's semi-empirical values for 64 Fe III 
lines from Ekberg (1993). Ekberg determined new wavelengths for the three lowest configurations 
of Fe III from laboratory measurements using a sliding spark discharge and spectrograph, with 
some energy levels computed assuming spin-orbit coupling, and transition probabilities calculated 
for all combinations. 

Data for the B III 25^8 - 2p'^P resonance doublet with lines at 2065.8 A and 2067.3 A, and 
hyperfine and isotopic components, are taken from Proffitt et al. (1999). The 2065.8 A, the 

stronger and less blended of the doublet, is our primary indicator of the boron abundance. This 
feature is blended with a weak Mn III line about 0.12 A from the ^^B line and 0.08 A from the 
^°B line. In our spectra, the B III and the Mn III lines are not resolved. Proffitt et al. 's high S/N 
and high resolution HST GHRS spectra of HD 35299 and HD886 clearly show the contributions 
of the ^^B, ^^B, and Mn III lines separately. Their synthesis results for these two stars suggest 
that the Mn III wavelength should be offset by 6 mA red ward of the Kurucz wavelength; also, 
they found that the abundance [Mn III] = —0.2 from that one feature, suggesting the gf-value is 
0.2 dex too large. We adopt the B III line list and transition probabilities given by Proffitt et al. 
(1999), and adjust the wavelength and oscillator strength of the blended Mn III line. [The weaker 
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B III line at 2067.3 A is blended with a strong Fe III line and weak Mn III line. The blending is of 
such severity that the 2067.3 A stellar line is not suitable for a boron abundance determination, 
although we note that the best boron fits to HD 216916, BD+56° 576, and the upper-limit for 
HD 36591 produce the best fit to the 2067.3 A B-Fe-Mn blend as weh (see Fig. g I, and D]. 

For all syntheses, we assume an isotopic ratio ^^B/^'^B = 4.0, the solar system ratio (Shima 
1963), consistent with the estimates given by Proffitt et al. (1999) from their line profile analyses of 
two sharp-lined B-type stars. We also consider a smaller ratio in our boron uncertainty estimates 
(discussed below). 

Proffitt & Quigley (2001) first noted the importance of the IS lines in this wavelength range, 
particularly one Cr II IS line that can come close to the B III A2065.8 feature depending on 
the stellar radial velocity. Since interstellar lines are quite sharp, we have identified their rest 
wavelengths from our spectrum of HD 205021. This star has a high S/N and rotation rate that 
make the broad stellar lines quite distinct from the sharp IS lines (e.g., see Fig. ^). We have 
identified four IS lines as Cr II (2055.60, 2061.58, and 2065.50) and Zn II (2062.01) from Morton 
(1991). We note the locations of the IS lines in our spectrum figures (see Fig.s |^ to where the 
IS lines appear at different wavelengths because the (stellar) rest frame is used. The IS line near 
the B III feature is only a problem in PG 0832+676. PG 0832+676 has two sets of IS lines, one of 
which appears to be blended with the Mn III-B III 2066 A feature. This affects the accuracy of 
our upper limit for boron in PG 0832+676, but nevertheless, the boron feature is clearly absent 
(see Fig. We also note in this figure that a Ti IV feature near 2067.6 A in PG 0832+676 is not 
well fit, but this is also seen in the spectrum of the hot star HD 34078, suggestive of atomic data 
uncertainties. 

Our final line list includes 2275 features between 2045 and 2078 A. All were included in the 
syntheses but many are negligible contributors. A few final fine adjustments were made to the 
line list. We examined the preliminary syntheses of three of the sharp-lined stars in this sample, 
HD 216916, HD 36591, and PG 0832+676, and made slight wavelength shifts to improve the final 
spectrum syntheses. These fine adjustments are reported in Table ^ We made no fine adjustments 
to the oscillator strengths. It is worth noting that our final line list does a remarkably good job 
at fitting the UV spectrum in our stars. This is noteworthy because UV line lists are notoriously 
incomplete and/or uncertain in their atomic data. Examination of our spectrum figures show very 
few missing lines and quite good fits, suggesting that the energy levels and transition probabilities 
are fairly accurate. 

3.2. Iron-group & Synthesis Parameters 

The iron-group abundances were determined from examination of a number of clean features 
in the spectra. Spectrum syntheses of these individual features were performed, tabulated, and 
the results averaged. This allowed us to minimize uncertainties due to poorly fit lines, as well as 



- 8 - 



lines which appear to have poor atomic data. In Table ^, we list the abundances relative to the 
meteoritic abundances from Grevesse & Sauval (1998), i.e., log(Fe) = 7.50 and log(Mn) = 5.53. 
The mean abundances in Table ^ were calculated by excluding line abundances that fall more 
than 2a (from the line-to-line scatter) from the mean. 

We estimate an uncertainty of approximately ±2 km s^^ in the macroturbulence, based on 
line profile fitting, but only ±1 km s~^ in microturbulence based on line strengths. Effects of these 
uncertainties on the iron-group abundances are quantified in Table |^ for a few representative stars. 
Uncertainties were found by fitting the iron-group features listed in Table |^ individually, then 
averaging the changes. Continuum placement and the adopted ^ and macroturbulence parameters 
dominate the errors in the iron- group abundances. NLTE effects are neglected throughout this 
iron-group analysis. Iron-group abundances are determined primarily from lines of the dominant 
species of the elements, e.g., Fe III, whose atomic statistical equilibrium rarely suffer from NLTE 
effects. 

We find the mean iron-group abundance of [Fe/H] = —0.06 ±0.15 from 6-23 features in seven 
main sequence B-type stars. This is similar to Cunha &: Lambert's (1994) results from 1-8 optical 
lines of Fe III in 16 Orion OBI B-type stars, [Fe/H] = —0.05 ±0.10, and also to Gies & Lambert's 
(1992) results from 1-3 optical hues of Fe II or Fe III in 31 Galactic B-type stars, [Fe/H] = ±0.22 
±0.20. Six stars have lUE spectra that were examined by Proffitt & Quigley (2001) for an overall 
scaling factor for heavy element abundances. Although they note that the scaling factor per star 
should not be interpreted as a precise determination of the heavy element abundance, we do note 
that their values are only about 0.2 dex lower than ours. 

Finally, due to the uncertainties in the B-type main sequence star temperature scale (discussed 
further in Section 4), we note that a slight reduction in the Gies &: Lambert (1992) temperatures 
(about 800 K, or 1 cj in the T^s uncertainty) will reduce our mean iron-group abundances by 
about 0.1 dex. 

3.3. Boron Abundances 

The LTE boron abundances listed in Table |^ are from spectrum syntheses, where we allowed 
the B III A2065.8 and Mn III A2065.9 line abundances to vary independently in order to achieve 
the best possible fit to the observation. Spectrum synthesis fits for all our program stars are 
shown in Fig.s |5|, and |6[ We did not attempt to constrain the Mn III A2065.9 line abundance, 
e.g., by setting the Mn abundance using other Mn III lines, because of uncertainties in the atomic 
data for other Mn III lines. We report the Mn III 2065.9 line abundance determined from the best 
synthesis of the B III blended feature in each star in Table ^ It is in good agreement with the 
iron-group abundances per star. 

To compute the uncertainty in the boron abundances, Alog(B/H), independent of uncertainties 
in Mn III 2065.9, we found it was necessary to fix the Mn III 2065.9 line abundance a priori. 
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Thus, we examined two methods (except for PG 0832+676, discussed below). Firstly, we simply 
set the value of Mn III from the best fit synthesis (in Table Q), then allowed only the B III line 
abundance to vary to calculate Alog(B). Secondly, we applied the mean iron-group corrections (in 
Table |^ to the Mn III line component, and allowed boron to vary to find a new best fit for each 
uncertainty parameter. We found Alog(B) to be very nearly identical from these two methods 
for all uncertainties investigated. The largest differences between the methods were <0.04 dex for 
A^Ma and A^. In Table ^, we report the uncertainties from the latter method. Only for AVrad 
did we vary B III 2065.8 and Mn III 2065.9 together (both line abundance changes are reported in 
Table |). 

For PG 0832+676, we could not fit the B III-Mn III blend directly because of an interstellar 
line occuring near the location of the Mn III line. For this star, we computed the iron-group 
abundances, and applied that result for the Mn III 2065.9 line abundance in order to determine a 
boron upper-limit. We used this method since the Mn III 2065.9 line abundance has been similar 
to the mean iron-group abundance in all of the other stars, i.e., [Mn/Fe] = —0.05 to —0.35, see 
Table Q for Mn and Table ^ for Fe. If we applied a lower Mn III 2065.9 line abundance (e.g., 
[Mn/Fe] = —0.3), then the change in the boron upper-limit is small (e.g., Alog(B/H) = +0.1). 

In summary. Table ^ shows that the most significant uncertainties in the boron abundances 
tend to be the continuum placement (thus, S/N of the data also) and the radial velocity and 
assumed ^^B/^'^B ratio (discussed further below) may be significant in some cases. Also, the Mn III 
2065.9 A atomic data could play an important role, but since the Mn abundance is not fixed for 
the boron syntheses, then uncertainties in the Mn atomic data will be compensated in the fit. 
Thus, the boron abundances should be accurate to ±0.1 dex. It is worth noting that lowering the 
temperature scale (discussed further in Section 4) has a negligible effect on the boron abundances 
(Alog(B/H) = 0.0 to -0.06). 

Only for the hottest star, HD 34078, is the derived boron abundance significantly sensitive 
to uncertainties in the atmospheric parameters. In fact, examination of the predicted equivalent 
widths for the B III 2065.8 A lines in Fig. ^ shows that the line strength is nearly constant between 
20,000 to 28,000 K, then drops dramatically at higher and lower temperatures. For example, near 
30,000 K, an uncertainty of ±2000 K can change the boron equivalent width by +50 %. From 
spectrum synthesis this corresponds to Alog(B/H) = +0.4, a substantial uncertainty for only a 
5 % uncertainty in T^s. For this reason, we examine our boron abundances versus temperature in 
Fig. |8|; we also include boron abundances from the literature on a homogeneous temperature scale 
(discussed below). There is no significant trend, but we mark the temperature limits where the 
B III 2065.8 A line should be the most reliable. 

Our LTE boron abundances in Table ^ are corrected for NLTE effects. For this paper, 
M. Lemke has extended his calculations (described in Cunha et al. 1997) to higher temperatures. 
The NLTE calculations were carried out for an abundance of 12 + log(B/H) =2.6. As seen in 
Fig. |9|, the NLTE corrections for the B III lines remain small, but reverse their sign in the hotter 
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stars. This is because the ground state of B III is overpopulated at the lower temperatures 
through colhsional couphng and ionization processes between B II and B III, and this is the only 
mechanism that produces B III NLTE effects at these temperatures. But as the temperature 
increases then the effects of B II on the ground state of B III are reduced. At the same time, the 
UV radiation field increases and there is a general shift towards B III as the main ionization stage 
such that photoionization dominates the B III statistical equilibrium. Eventually, this causes the 
sign reversal in the B III NLTE corrections. 

Boron is clearly present and undepleted in only two stars in our sample: BD+56° 576 and 
HD 216916. That the boron abundance is quite different in stars of very similar atmospheric 
parameters is obvious from inspection of the spectra of HD 216916 and HD 16582. For example, 
compare the large differences in the B III-Mn III 2066A blend in Fig. ^ for these two stars with 
insignificant differences in their other features (also see Fig.s[l|, |2|, and^. Also, both of these stars 
are (5 Cep variables (see Table ^) , thus we conclude that the (3 Cep variability mechanism is not a 
cause of surface mixing. Similarly, the Orion star HD 36591 has similar atmospheric parameters 
to another Orion star HD 35299, yet vastly different boron line strengths are seen in Fig. |^. 

Six of our stars have lUE spectra that were examined by Proffitt & Quigley (2001). Owing to 
the lower S/N and resolution of lUE spectra, their boron abundances are less precise than ours. 
When boron is strong, as in HD 216916, our results are in excellent agreement and with similar 
uncertainties (±0.2 dex). On the other hand, Proffitt & Quigley report detecting the boron line in 
HD 34078, which we do not confirm (see Fig. P). For the other four stars in common, we report 
significantly lower abundances (HD 16582) or lower upper-limits. 

4. Discussion 

4.1. The Initial Abundance of Boron 

In order to pursue our principal goal of testing predictions of boron-depletion due to 
rotationally induced mixing, it is necessary to establish the initial boron abundance of these 
local B-type stars. (PG 0832+676 is discussed separately.) We suppose that the highest boron 
abundances across our sample and others provide a close approximation to the initial boron 
abundance. We further assume that stellar boron abundances that are substantially less than the 
presumed initial abundance do not reflect an unusually low initial value, but rather arise from 
internal processes that deplete boron. 

The mean NLTE boron abundance from BD+56° 576 and HD 216916 is 12+log(B/H) = 2.3. 
Proffitt et al. (1999) 's analysis of another two stars gives 12+log(B/H) = 2.4 (considering NLTE 
corrections), and the mean NLTE abundance of stars with boron detections from the sample 
considered by Proffitt & Quigley (2001) gives 12+log(B/H) = 2.4. Cunha et al. (1997) found near 
solar boron abundances from two B-type stars in the Orion OBIc subassociation, although those 



- 11 - 



results were from the B II 1362 A line. Boron from the B I 2497 A line in F-G main sequence stars, 
with -0.15 > [Fe/H] < +0.15, is 12+log(B/H) = 2.6 from Cunha et al. (2000a). 

These stellar abundances should be closely related to the interstellar abundance of boron, as 
standard evolutionary models, as well as those including rotational effects (see below), predict 
survival of surface boron through the pre-main sequence phase. Boron is detected in interstellar 
diffuse clouds by the B II 1362 A resonance line (Jura et al. 1996; Lambert et al. 1998; Howk, 
Sembach, & Savage 2000). The lower limit to the interstellar abundance, 12+log(B/H) > 2.4 
suggested by Howk et al., is consistent with the above stellar abundances. The lower limit is given 
in recognition of the possibility that boron may be depleted onto interstellar grains. It is most 
unlikely that boron is present in grains, if, as is widely supposed, grains are formed primarily in 
stellar outflows from red giants whose outer envelopes will be highly depleted in boron. 

These stellar and nebular abundances are slightly lower than the solar system (meteoritic) 
value of 12+log(B/H) = 2.78 (Zhai & Shaw 1994). For the B-type stars and the interstellar gas, 
this reduction is in line with similar underabundances (relative to the Sun) for oxygen and other 
light elements. Since our primary interest lies in the B-type stars depleted in boron by an order 
of magnitude or more, the difference between an initial boron abundance of 2.4 and 2.8 is largely 
irrelevant. 

It is also worth investigating the boron abundance in the Orion association as a fiducial 
point, since there are both stellar and interstellar determinations available. In general, 
2.4 < 12+log(B/H) < 2.9 from B-type stars and interstellar sight lines (among the references 
above). Some stars in Orion show large boron depletions though. Two of our stars are or were 
members of the Orion association: HD 36591 belongs to the lb subgroup, and HD 34078 (AE Aur) 
is a well known runaway star. Neither star shows a detectable B III line. While the boron depletion 
in HD 34078 might be related to its transformation as a runaway star (e.g., contamination of 
the star by ejecta from its exploding companion), HD36591's low boron abundance is plausibly 
identified as arising from internal effects (discussed below). This result is significantly different 
from other boron determinations from B-type stars in Orion (Cunha et al. 1997, Proffitt et al. 
1999, Lemke et al. 2000). 

Low boron abundances have been found in some other stars in Orion. Cunha et al. (1999, 
2000b) determined the boron abundances in three G-type stars in Orion from the B I 2497 A line. 
One of these, BD— 5° 1317, resulted in an abundance 5x less than solar, 12+log(B/H) = 2.1 ±0.2. 
Surprisingly, this star also proved to be oxygen rich (Cunha et al. 1998), 3x larger than the Orion 
nebular abundance (Esteban et al. 1998). These results lead Cunha et al. (2000b) to propose 
a boron-oxygen anticorrelation, possibly the result of dilution of interstellar boron by relatively 
boron-poor but oxygen-rich ejecta from recent Type II SNe. 

Proffitt & Quigley (2001 = PQOl) also discovered severely boron-depleted B-type stars, which 
muddles the issue of a boron-oxygen anticorrelation, particularly as the boron depletions are not 
related to oxygen enrichment. In Fig. |l^ (left panel), we show the boron-oxygen anticorrelation 
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from published data for Orion, as well as our star HD 36591. (We note that PQOl include seven 
more Orion stars with depleted boron, however five are hot stars like HD 34078, which we consider 
less reliable boron abundance indicators, and two have no published oxygen abundances). In the 
right panel, we include only stars where boron is determined from the reliable B III feature (the 
boron and oxygen abundances plotted are from Table ^, which is discussed further below). The 
existence of a boron-oxygen anticorrelation is not obvious from the B III data alone, although we 
note that much of the B III data is from lower quality lUE spectra. 

The phenomenon of boron-depletion is not confined to the Orion association. Proffitt 
et al. (1999) found that the field star HD 3360 is boron-depleted, and we have found several 
boron-depleted stars in other OB associations, as did PQOl. Fliegner, Langer, & Venn (1996) 
first proposed that boron depletions may result from rotationally-induced mixing during the main 
sequence lifetime of B-type stars. We discuss this in the next section. 

4.2. Rotationally-induced Mixing 

Recent models by Heger et al. (2000) follow the evolution of the angular momentum 
distribution and the occurence of associated mixing processes in massive stars from the pre-main 
sequence through core collapse. They find that rotational mixing can affect stellar surface 
abundances, stellar life times, and the evolution of a star across the HR diagram (Heger &: Langer 
2000 = HLOO). These new models have fundamental consequences for the age determination of 
young stellar clusters, and they can reconcile the long-standing mass discrepancy problem (that 
spectroscopic masses differ from those inferred from stellar evolution tracks, Herrero et al. 1992). 

In general, a rotating star has a lower effective gravity, thus it acts like it has less mass at 
core-H ignition. Later, during core-H burning, rotationally induced mixing of protons from the 
envelope into the convective core and of helium from the core into the envelope will lead to higher 
luminosities compared to non-rotating models. Also, the effective temperature is lower, thus the 
ZAMS position on the HRD changes from a non-rotating model. Evolution on the main sequence 
then depends on the star's mass and rotation rate, as well as the efficiency of mixing in the upper 
convective core and stellar interior. 

Heger & Langer 's models cover a mass range of 8 to 25 Mq, and a range of rotational 
velocities from to ~450 km s~^. They computed two sets of models, with different assumptions 
on the efficiency of rotational mixing in layers containing a gradient in the mean molecular 
weight n (called a /i-barrier). In one set, the ^-barriers are ignored, and the overall efficiency of 
rotational mixing had to be reduced in order to meet observational constraints. In the other set, 
/i-barriers are taken into account, and again the overall mixing efficiency above the /i-barrier is set 
by observational constraints (see Heger et al. 2001 for details). The differences between the two 
sets of models reflects the remaining uncertainties in the theoretical description of the rotational 
mixing processes. Finally, it is interesting to note that even though the rotationally induced 
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/i-barrier inhibits mixing just above the core, the current models show that the envelope above the 
/i-barrier is very well mixed, thus CN-processed gas is more evident up through the photosphere 
then when the ^-barrier is ignored. 

The range in velocities examined by the models is thought to be typical for B-type stars. For 
example, Fukuda's (1982) statistical study of rotational velocities found B2 V and B2 IV stars 
with mean vsini values of 154 and 110 km s~^, respectively, and that OB stars generally span 
a range in vsini of 100 to 400 km s^^. Also, de Jager (1980, p. 50) summarizes mean equatorial 
velocities for B1/B2 V stars (the majority of stars in this paper) as ~200 km s~^. Although all 
the stars in this paper, and others with good abundance analyses, are sharp lined objects, the 
usual assumption is that many of these stars are (near) pole-on rotators. 



In Fig. 11, we show the HLOO models for and 200 km s~^, with and without ^u-barrier effects 
in the log g-logTefj diagram. We also include our stellar targets, and more B-type stars with 
boron and CNO abundances from the literature, to show their approximate mass range and main 
sequence ages (differences between the symbols are discussed below). We notice that the model 
ZAMS is offset from the highest gravity stars. This is most likely an artefact of uncertainties in 
the temperature scale (discussed below), which subsequently affects the determinations of gravity 
and ^. 



4-2.1. Boron versus Nitrogen: Predictions 

Heger & Langer's (2000) models can reproduce previously unexplained observations of massive 
stars (e.g., various abundance anomalies) and the models make certain predictions (see their 
discussion of the observational evidence and testable predictions). The prediction that concerns 
us in this paper is the relationship between boron and nitrogen^ during main sequence evolution. 

Mechanisms that can deplete boron at the surface of stars can also enrich the surface in the 
products of the CN-cycle; material that is rich in nitrogen and depleted of carbon, yet where the 
sum C-|-N is preserved. Extremely deep mixing could also bring ON-cycled material to the surface, 
preserving the sum C-|-N-|-0. Thus, in Fig. we show the boron-to-nitrogen relationship and the 
rate of boron depletion predicted by various HLOO models for main sequence evolution (through 
core-H burning). We show models that include /i-barrier effects for 12 Mq at five rotational 
velocities, and for comparison we have included two 15 models at 200 km s~^ and 450 km s~^ 
as well as the 10 Mq and 15 Mq models without the //-barrier effects. These figures allow us to 
examine the effects due to rotation rate, mass, age, and mixing efficiencies. 

One of the most interesting results is that the boron-nitrogen relationship is nearly the same 
for all models. This occurs because boron is destroyed throughout most of the star very early in 



* We note that boron versus the ratio N/C is also a good indicator, but our opinion is that the small depletions in 
carbon due to mixing do not clearly compensate for the added uncertainties in the B-type star carbon abundances. 
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its main sequence lifetime (<10^ yrs). Only the outermost layers of the stellar interior are cool 
enough to retain the initial boron (the outer ~1 Mq for 10-20 Mq stars). Meanwhile during this 
early ZAMS phase, the CN-cycle begins in the core, transforms carbon into nitrogen, and mixes 
this gas into the interior through convective overshoot and/or semi-convection. As the star ages 
on the main sequence, rotation will mix the CN-cycled gas up through the stellar interior, while 
also mixing the pristine outer layers downwards. Fig. |l^ illustrates that mixing the radiative 
envelope of a massive main sequence star will always have the same consequences for the surface 
abundances; boron is depleted first, and nitrogen is enriched only later (also shown in Fliegner 
et al. 1996). Due to the similarity in the structure of the envelopes of the stars considered, in 
particular the shallowness of the surface layer that contains pristine boron and the depth of the 
layers were nitrogen is enriched, then the same abundance changes are predicted independent of 
all parameters. 

On the other hand, the timescales for the B-N abundance changes do differ significantly 
between the models. Since rotational mixing is less effective for slower rotation rates or lower 
masses, then the abundance changes simply take longer to occur. Furthermore, very few models 
predict high N/C ratios at the end of the main sequence lifetime, thus additional processing/mixing 
after the main sequence phase (such as the first dredge-up) can be examined through the N/C 
ratios of more evolved stars. 

We comment once again that mass loss could produce the same B-N abundance changes as 
rotation (having the same ZAMS starting conditions as well). However, the mass loss rates for 
B-type stars are much smaller than required (as discussed earlier), thus favoring rotation for this 
phenomenon. Also, boron depletion and nitrogen enrichment are predicted by mass transfer in 
a close binary system. Wellstein (2000) and Wellstein et al. (2001) predict changes in surface 
abundances of boron and CNO due to the transfer of nuclear processed matter. The current models 
predict much larger boron depletions though, and there is always a very strong CNO-signature as 
well. Therefore, the existence of moderately boron-depleted stars without nitrogen enrichments is 
a unique signature of rotational mixing effects. 

4-2.2. Boron versus Nitrogen: Observations 

In the abundance analysis described above, we have adopted the Gies & Lambert (1992 
= GL92) atmospheric parameters. Systematic differences between GL92 and other B-type star 
analyses suggest that this temperature scale is too hot (e.g., Korotin et al. 1999, Cunha & 
Lambert 1994 = CL94). In fact, GL92 did increase their temperature scale by 3.4% from the 
photometric determinations. This increase removed an unexpected relationship between NLTE 
CNO abundances and temperature. However, this relationship appears to have been due to the 
NLTE CNO abundances having been determined from Gold (1984) model atmospheres, instead of 
the more heavily line-blanketed Kurucz models. 
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In Table we have reduced the GL92 temperature scale by 3.4%, and applied their NLTE 
CNO abundance corrections listed in their Table 9 (the A values). We have also applied a 
correction to account for the Gold-Kurucz offsets, as tabulated by CL94 (their Table 10). (Note 
that this latter correction is derived from differences in LTE abundances between Gold and Kurucz 
models, and we assume the same will apply to the NLTE abundances, as assumption that should 
be checked). The GL92 temperatures are now in good agreement with CL94 values for stars in 
common. In the case of BD+56° 576, which was not analysed by GL92 or CL94, we take the 
parameters and NLTE CNO results of Vrancken et al. (2000), but add 0.4 dex to their published 
nitrogen value to account for the systematic differences they discuss between their analysis and 
those by GL92 and Kilian (1992, 1994). 

The GL92 NLTE nitrogen abundances are now in excellent agreement with CL94 values for 
stars in common. Only the hottest star, HD 36960 shows a non-negligible nitrogen difference of 
~0.1 dex. However, it is surprising that the carbon and oxygen abundances often differ between 
these analyses, by up to 0.25 dex. For carbon, the abundances differ when GL92 include two 
strong C III lines. For oxygen, the differences are traceable to the adopted ^ values. CL94 adopt 
higher ^ values as their analysis is more sensitive to ^ due to some stronger lines. For consistency, 
we adopt the corrected GL92 values for CNO in Table ^ whenever possible. 

The boron abundances determined from the B III feature are not sensitive to these 
uncertainties in temperature, except for the hottest and coolest stars that we consider less reliable 
boron indicators (discussed above), and thus that we have excluded from Table ^ (including 
our hot target HD 34078). We also note that the uncertainties in Teff and/or gravity do not 
significantly affect the boron-nitrogen relationship that we seek; i.e., the B III and N II line 
abundances react similarly to the atmospheric parameters in most cases (only in the coolest stars 
is nitrogen more sensitive to the atmospheric parameters). Thus, uncertainties in the parameters 
cannot induce a B-N anticorrelation. 



4-2.3. Predictions versus Observations 

To compare the boron and nitrogen abundances in B-type stars to the model predictions, 
we have gathered abundances from this analysis and the literature, corrected them to be on 
a homogeneous temperature scale, and applied the known Gold-Kurucz model atmosphere 
abundance corrections (discussed above). We have also limited the atmospheric parameter range 
of our sample to 18,000 K < Teg < 29,000 K, where the B III line strength is in a plateau, and 
log g >3.4 to examine main sequence stars only. GL92 showed that their sample of B-type stars 
contained some stars mildly depleted of carbon and enriched in nitrogen with no measureable 
spread in the oxygen abundances. The C-N signature is indicative of CN-cycled material in the 
atmospheres of some stars. There is no evidence for ON-cycled gas in main sequence B-type stars. 



In Fig. 13, we can see that most stars gather around the point 12-|-log(B/H) = 2.6 and 
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12+log(N/H) = 7.7, consistent with the ISM boron abundance in Orion (discussed above) and the 
ISM nitrogen abundance (e.g., 12+log(N/H) = 7.8, Esteban et al. 1998). Possibly, there is a small 
intrinsic spread - initial boron may be higher (about 2.8) and initial nitrogen lower (about 7.6) - 
or perhaps some small depletion by rotation has already occured. Boron and nitrogen predictions 
from HLOO for the 12 Mq models (with 200 km s~^ and //-barrier effects) are shown in Fig. [l^. 
The same model is shown scaled with three different initial boron and nitrogen abundances (all 
initial abundances used here are within the observational uncertainties) to better compare to the 
stellar observations. Only the predicted ^^B abundances are traced. The isotope ^'^B is more 
readily destroyed by protons, such that one would assume that the "B/^°B would increase m 
regions of partial boron destruction. However, the HLOO models predict a slower depletion rate for 
^'^B. Thus, ^^B/^'^B = 4 adopted for the ZAMS is predicted to reduce to >2 in the stellar envelope 
by the end of the main sequence. We note that these differences in the boron isotopic ratio have 
a negligible effect on the boron abundance determinations in most stars (see Table ^). Only two 
stars have low-boron determinations, HD 16582 (this paper) and HD3360 (Proffitt et al. 1999). 
Several others have interesting upper limits that follow the boron-nitrogen trend predicted by 
the rotating models. HD 36591 in Orion is particularly interesting since it shows a strong boron 
depletion with no enrichment in nitrogen; the model predictions are most consistent with this star 
for the lowest initial nitrogen (=7.6) abundance. 



In Fig. |14|, we examine the boron abundances with stellar age. Three stars in this paper 
are in young clusters, and show boron depletions that are fit well by rotating model predictions. 
Two of these stars have nitrogen enrichments, thus they could also be explained by binary 
accretion (although only one is in a known binary system), but the third has unenriched nitrogen. 
As previously stated, the best interpretation for boron-depletion without nitrogen-enrichment 
is rotational mixing. This third star is HD 36591, again, and we find that the most extreme 
models, e.g., 20 Mq, fe(j=450 km s~^, reproduce its age, boron, and nitrogen abundances well. 
Unfortunately, the atmospheric parameters for this star make it very unlikely that it is a 20 M© 



star (see Fig. 11). Perhaps mixing efficiencies above /x-barriers can be higher then currently 
predicted. 



Also seen in Fig. 14 is that some B-type stars in older clusters have normal boron abundances. 
There are two possiblities for this. First, they may be true slow rotators (not rapidly rotating 
stars seen pole-on as is usually assumed for sharp-lined objects). This is likely true for some 
stars, e.g., HD 216916 is an eclipsing binary (Pigulski & Jerzykiewicz 1988), thus its low vsmi 
is probably close to its equatorial velocity. Also, some older subgroups of OB associations have 
been reported to have an excess of slow rotators (Guthrie 1984, also see the discussion by Grebel 
et al. 1992). Second, some stars may be lower mass objects, e.g., the 10 Mq models require more 
time to mix hotter gas from the interior to the surface, thus nitrogen and boron remain near 
their initial abundances longer. Schoenberner and Harmanec (1995) determined masses of main 
sequence B-type stars in binary systems, and found that stars in our primary temperature range 
(24000-28000 K) have masses of 10-12 M©, consistent with Fig. |l|. 



-17- 



In Fig. 11, we also examine the boron (from B III only) and nitrogen abundances on a 
Teflf-gravity diagram to examine the abundance distributions across the main sequence. The stellar 
data is divided into three groups: normal stars (B>2.2 and N<7.8), boron depleted stars (B<2.2 
and N<7.8), and nitrogen rich stars (B<2.2 and N>7.8). As expected, there are no N-rich and 
B-normal stars. Evolution tracks from HLOO are shown for 8-20 M0 through the core-H burning 
phase. In general, the three groups of stars have the same range in T^s and gravity, and thus in 
age and mass, suggesting that their differences are due to another parameter, like rotation rate. 
It is interesting that most of the B-depleted/N- normal stars are near the ZAMS, as is predicted 
for rapid rotators. Also, one might argue that among the most massive stars, more of them show 
B-depletions; this would be consistent with the models if the rotation rates are similar throughout. 

Finally, it is pleasing that BD+56° 576, the star with the lowest surface gravity of our sample, 
has apparently retained its initial complement of boron and has unenriched nitrogen. This star is 
in an eclipsing binary system, which implies its t;sini value is very close to its current rotational 
velocity. Thus BD+56° 576 is probably a slow rotator today, and was probably a slow rotator on 
the ZAMS (although it may have spun down somewhat from its ZAMS rotation rate). Only one 
other star shows unenriched nitrogen and a lower surface gravity, HD 30836, but PQOl show that 
this star has depleted boron, suggestive of rotational mixing. (This star is in Ori OBI, but it is 
not in Fig. 14 since its subassociation, thus age, is not clear). 



4.3. PG 0832+676 

Our determination of the metalhcity of PG 0832+676 is [Fe/H] = -0.88 ±0.10, which 
compares with [Fe/H] = —0.51 given by Hambly et al. (1996) from a single optical Feiii line in 
a differential analysis with respect to the Galactic B-type star HR 1886. Combining our [Fe/H] 
with the mean a-element (Mg, Si, and S) abundance from Hambly et al. , we obtain [a/Fe] = 0.5, 
which is the expected value for a star of this metalhcity to within the errors of measurement (e.g., 
[0/Fe] = 0.3 for [Fe/H] = —0.9, as for the mildly metal-poor clusters M4 and M71, Ivans et al. 
1999, Sneden et al. 1994, although the 0/Fe ratio in metal-poor stars is currently controversial, 
c.f., Lambert 2000). 

Relative to [Fe/H] = —0.88, PG 0832+676 is enriched in carbon, nitrogen, and oxygen, 
according to Hambly et al.'s results: [C/Fe] ~ [N/Fe] ~ [0/Fe] ~ 0.7 to 0.8. This 0/Fe ratio is 
somewhat larger than expected for a normal metal-poor star in the halo (see a/Fe ratio comments 
above), and carbon and nitrogen appear enriched relative to normal stars. These results are 
consistent with the identification of PG 0832+676 as a post-AGB star with carbon and oxygen 
added in the AGB phase, and nitrogen enhanced by the earlier first dredge-up. Certainly this 
scenario is consistent with our non-detection of boron. Our limit of 12+log(B/H) < 0.60 is much 
less than the expected initial abundance of 12+log(B/H)=1.9 at [Fe/H] = —0.9 (based on the 
analysis of 14 cool dwarfs by Cunha et al. 2000a). 
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Thus, in summary, the chemical pattern in PG 0832+676, low [Fe/H], high [a/Fe] and 
[CNO/Fe], very low [B/H], suggests that this star is a post-AGB star, in agreement with Hambly 
et al. 's (1996) conclusion that PG 0832+676 is a highly evolved star. 



5. Conclusions 

Boron in hot stars, like lithium in cool stars, is shown to be a tracer of some of the various 
processes affecting a star's surface composition that are not included in the standard models of 
stellar evolution. If the initial boron abundances of local hot stars are similar from star-to-star and 
association-to-association, then the large spread in boron abundances, a factor of at least 30 across 
our sample, shows that boron abundances are a clue to unravelling the non-standard processes 
that affect young hot stars. In this paper, we have focussed on the role of rotationally-induced 
mixing. 

Models of stars with masses near 10 Mq (Heger & Langer 2000) show that rotationally-induced 
mixing during main sequence evolution can reduce the surface boron abundance and increase 
surface nitrogen. A signature of rotationally-induced mixing is that the initial decline of boron 
precedes an observable change in nitrogen. The correlation between the decline of boron and the 
rise in nitrogen is almost independent of rotational velocity and mass of the star. On the other 
hand, the extent and rate of change of the abundances are dependent on both velocity and mass 
(Fig- 0)- 

Our results confirm Proffitt et al. 's (1999) discovery that boron may be quite severely 
depleted in otherwise normal B-type stars; 6 of the 8 stars in our program show a reduced boron 
abundance. That boron depletions are not a rare occurrence is consistent with Proffitt & Quigley's 
(2001) survey of lUE spectra. Boron depletions are attributed to rotationally-induced mixing 
(although N-rich stars could also be explained by binary mass transfer). The correlation between 



the boron and nitrogen abundances follows the predicted trend quite well (Fig. 13). 



Our use of stars from various OB associations allows us to investigate the rate of change of 
the surface abundances. Stars showing a normal boron abundance are found at all the investigated 
ages; this is consistent with predictions for rotationally-induced mixing for low equatorial velocities 
(100 - 200 km s-i). The observed vsini velocities are considerably smaller than these limits, but 
we assume most stars are rapidly rotating yet seen pole-on. Two stars appear to have depleted 
boron at a far faster rate than expected; HD 36591 and HD 205021 are very young stars in our 
program with no detectable boron. In Fig. they appear close to the locus for stars of 20 



Mq rotating at 450 km s~^, but their locations in Fig. 11 imply that the stars are very early 
main sequence stars with masses of about 12 — 13 Mq. This discrepancy may indicate that 
rotationally-induced mixing has been underestimated at high rotational velocities or that it is 
more efficient then expected at lower masses. 



The spectral window around the B III resonance lines proves well suited for a determination 
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of the abundance of the iron-group elements. For the local B-type stars in our sample, we show 
that they have an approximately solar metallicity. For the halo B-type star PG 0832+676, we find 
a low metallicity and an absence of boron which suggest the star is probably a post-AGB star. 

Further insights into rotationally-induced mixing will require STIS spectra of additional 
rapidly rotating stars. A few stars examined by Proffitt & Quigley (2001) have low boron 
upper-limits and unenriched nitrogen, and are ideal targets for follow-up studies. Additional stars 
that are sharp- lined (low vsini), but rapidly rotating, are also suitable; however, in the absence 
of a way to extract the angle of inclination, information about abundances and rotation will 
require a statistical treatment. We also note that stars hotter than the present sample are not 
appropriate; the B III lines become too weak (see Fig. 0) and boron, present as the He-like ion, is 
lost to spectroscopic scrutiny. 
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Table 1. HST STIS Observing Information for Galactic B-stars 



btar 


V 


Gr at /Slit 


Exposure(s) 


T— \ 1 

Date 


o /at 

b/N 


BD+56 576 


9.38 


E230M 


1680s at 


Ac2124 


11 FEB 99 


100 






0.2x0.2 


+1280S at 


Acl978 












+1509S at 


Ac2269 






HD 16582 


4.07 


E230H 


1634s at 


Ac2113 


23 JAN 99 


90 






0.2x0.05ND 


+1140S at 


Ac2063 












+1140s at 


Ac2063 












+1368S at 


Ac2013 












+1328S at 


Ac2013 






HD 34078 


6.00 


E230H 


432s at 


Ac2063 


15 MAR 00 


50 






1x0 03 


+432s at 


Ac2013 






HD 36591 


5.34 


E230M 


432s at 


Ac2124 


09 FEB 99 


100 






0.2x0.05ND 


+432s at 


Ac2124 












+432s at 


Ac2124 












+432s at 


Acl978 












+432s at 


Acl978 












+432s at 


Ac2269 












+432s at 


Ac2269 






HD 50707 


4.83 


E230M 


216s at 


Ac2124 


06 NOV 98 

\J\J -1. 1 V_x V t_/ K-J 


75 






2x0 05ND 


+216s at 


Acl978 












+216s at 


Ac2269 






HD 205021 


3.23 


E230H 


900s at 


Ac2063 


19 FEB 99 


100 






2x0 05ND 


+900s at 


Ac2063 












+684s at 


Ac2013 












+684s at 


Ac2113 












+684s at 


Acl963 






HD 216916 


5.59 


E230M 


1260s at 


Ac2124 


16 DEC 98 


130 






0.2x0.05ND 


+262s at 


Ac2124 












+338s at 


Ac2124 












+888s at 


Acl978 












+972s at 


Ac2269 






PG 0832+676 


14.15 


E230M 


3360s at 


Ac2124 


08 DEC 98 


100 






0.2x0.2 


+3360S at 


Ac2124 


cvz 










+3360S at 


Ac2124 












+3300S at 


Ac2124 
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Table 2. Variability and Binary Information from the Literature 



Star 


Name 


Variable 


Binary 


Gluster 


CI. Age 
(Myr) 


REFS 


BD+56 576 






eclipse 


X Per 


11.5 


KP97, 037, T+84 


HD 16582 


5 Get 


(3 Gep 




Gas-Tau 


~50 


H+94, G63, deZ+99 


HD 34078 


AE Aur 


var 


visual? 


Ori OBI 




G87 (run-away) 


HD 36591 




var 


visual 


Ori OBIb 


1.7 ±1.1 


077, BV97, WH78, B+94 


HD 50707 


15 GMa 


P Gep 




Goll 121 


~5 


H+94, deZ+99 


HD 205021 


/3 Gep 


P Gep 


spect 


Gep OBI? 


2 ±1 


T+97, F69, M+95 


HD 216916 


16 Lac 


P Gep 


eclipse 


Lac OBI 


12-16 


PJ88, JP99, F69, B91, deZ+99 


PG 0832+676 








Halo 







Note. — References are B91 = Blaauw 1991, BV97 = Brown &; Verschueren 1997, B94 = Brown et al. 
1994, G63 = Grawford 1963, deZ+99 = de Zeeuw et al. 1999, F69 = Fitch 1969, G87 = Gies 1987, H+94 
= Heynderickx et al. 1994, JP99 = Jerzykiewicz & Pigulski 1999, KP97 = Krzesihski & Pigulski 1997, 
M+95 = Massey et al. 1995, 077 = Olsen 1977, 037 = Oosterhoff 1937, PJ88 = Pigulski & Jerzykiewicz 
1988, T+84 = Tapia et al. 1984, T+97 = Telting et al. 1997, WH78= Warren & Hesser 1978. 
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Table 3. Atmospheric Parameters from the Literature 



Star 


SpType 




log g 


f sini 




REF 






(K) 




(km s~^) 


(km s~^) 














(NLTE) 




BD+56 576 


B2III 


22500 


3.40 


<50^ 


9 


v±oo 




B2III 


21500 


3.6 




12 


L88 


HD 16582 


B2IV 


23750 


4.08 


15 


3.8 ±5.1 


GL92 


HD 34078 


09.5V 


31420 


4.07 


27 


8.0 ±0.6 


GL92 


HD 36591 


BlIV 


27380 


4.15 


11 


4.4 ±3.7 


GL92 




BIV 


26330 


4.21 




9 


CL94 


HD 50707 


BlIV 


27710 


4.04 


39 


4.2 ±4.6 


GL92 


HD 205021 


BlIV 


26740 


4.16 


28 


3.3 ±3.7 


GL92 


HD 216916 


B2IV 


24050 


3.90 


13 


3.6 ±4.6 


GL92 


PG 0832+676 


BIV 


23000 


3.7 


1 


yb 


H96 



Note. — References are GL92 = Gics & Lambert (1992), V+OO = Vrancken 
et al. (2000), H96 = Hambly et al. (1996), CL94 = Cunha & Lambert (1994), 
L88 = Lennon et al. (1988). For BD±56 576 and HD 36591, the parameters 
from V±00 and GL92 were adopted, respectively. 

'^For BD±56 576, vsini = 17 km s^^ from our spectrum syntheses, 
in PG 0832+676 from a LTE analysis. 
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Table 4. Boron and Iron-Group Abundances from STIS Spectroscopy 



Star 


RV 






[M/H] log(5/H) log(S/H) log(Mn III) 
LTE NLTE A2065.9 


BD+56 576 


1 


20 


4 


-0.16 ±0.17 


2.43 


2.25 


5.24 


HD 16582 


12 


12 


2 


-0.15 ±0.15 


1.27 


1.16 


5.31 


HD 34078 


56 


7 


3 


-0.2 ±0.2 


<2.0 


<2.2 


5.1 


HD 36591 


31 


18 


2 


±0.02 ±0.12 


<1.45 


<1.38 


5.43 


HD 50707 


28 


20 


6 


±0.1 ±0.1 


<1.5 


<1.5 


5.4 


HD 205021 


-29 


12 


2 


-0.16 ±0.23 


<1.00 


<0.90 


5.00 


HD 216916 


-18 


19 


1 


±0.16 ±0.18 


2.44 


2.31 


5.33 


PG 0832+676 


-68 


20 


2 


-0.88 ±0.10 


<0.75 


<0.60 


4.65 



Note. — Abundances have been determined from spectrum syntheses using the 
model atmosphere parameters listed here and in Table ^ Radial velocities, ^, and E,Ma 
(macroturbulence) values are determined from the iron-group features. We estimate ARV 
and IS.^Ma ~2 km s~^ based on line profile shapes, and <1 km s~^. The B III 
A2065.8 and Mn III A2065.9 abundances were allowed to vary independently for the best- 
fit syntheses, thus we report the Mn III 2065.9 line abundance here. 
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Table 5. Iron- Group Wavelength Offsets 



iLiiem 


\ Zt/'t T"r> \ 

A(KUKj 


A(JNll,Wj 


Mn III 


2048 949 


2048 909 


Mn III 


2049 357 


2049.314 


Mn III 


2049.682 


2049 663 


Mn III 


2052 739 


2052 673 


Mn III 


2063 337 


2063 397 


Mn III 


2065 886 


2065 892* 


Fe III 


2050.743 


2050.738 


Fe III 


2052.271 


2052.250 


Fe III 


2053.524 


2053 508 


Fe III 


2054.492 


2054.475 


Fe III 


2055.863 


2055.848 


Fe III 


2056.152 


2056.147 


Felll 


2057.059 


2057.052 


Felll 


2058.209 


2058.199 


Felll 


2058.566 


2058.558 


Fe III 


2059.692 


2059.678 


Fe III 


2064.980 


2065.030 


Fe III 


2068.896 


2068.983 


Felll 


2070.976 


2070.976 


Fe III 


2076.322 


2076.309 


Ni III 


2045.412 


2045.387 



Note. — Our line list originates from Kurucz 
(1998, CD18), with wavelengths for eight Fe III 
lines updated from Proffitt et al. (1999) and 64 
others updated from Ekberg (1993; wavelengths 
and (^/-values). In this tabic, we note additional 
adjustments to some line wavelengths based on our 
preliminary syntheses of sharp-lined stars in our 
sample. One line from Ekberg, 2070.976 A, was 
moved back to the original Kurucz wavelength, but 
the new gf-value retained. 

*New wavelength from Proffitt et al. (1999), who 
also suggest log5/=— 0.241. This feature is blended 
with the B III 2065.8 line. 
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Table 6. Iron-group Abundance Results 



A 


blement(sj 


[M/HJ 


[M/HJ 


r 71 /T /TTl 

[M/HJ 


[ 1\ /T /TTl 

[M/HJ 


[ i\ /r /TTl 

[M/HJ 


r 71 /T /TTl 

[M/HJ 


\ 1\ /T /TTl 

[M/HJ 


r 71 /r /TTl 

[M/HJ 


(A) 




BD56 


16582 


34078 


36591 


50707 


205021 


216916 


PG0832 


904^ 3Q 


Ni TTT 




—0 07 




-1-0 OQ 






-1-0 16 




9048 Q1 


Mn TTT 

1VJ.11 111 


00 


-1-0 03 




00 


-1-0 9^ 


—0 05 


+[) 08 


—0 Qfi 


9n4Q 37 


Fp TTT-I-A/Tn TTT 

± C llin^-LV±ll 111 


00 


—0 94 


—0 3^ 


-1-0 05 


-1-0 9^ 


-1-0 95s 


-1-0 90 


—0 QO^ 


9n4Q fifi 


A/Tn TTT 

IVlll 111 




-/-/7 


—0 3^ 


-1-0 05 


-1-0 9^ 


-1-0 95s 


-1-0 77 


-1-0 OQ 




T7p tTT-I-Fp TV 
± c iiin^i c- 1 V 


—0 90 


—0 30 




-1-0 95 










90^9 9^ 


Fp TTT 

A. C 111 


—0 90 


—0 95 




-1-0 1 3 






-1-0 35 


—0 Q3 


i*»JtJ*J . (J J- 


Fe III 


— fl 5R 


—0 35 


—0 4 






00 


-1-0 1 


—0 98 


9054 5fi 


Fp TTTx3 


—0.24 


—0 06 


—0 4 


00 





-1-0 95 


-1-0 1 q 


—0 84 


9055 85 


Fe III 




-4-0 9^ 


—0 2^ 


-1-0 08 


+0 9^ 






— OP/, 


905fi 1 5 


Fp TTT 




-1-0 09 


—0 2^^ 


-1-0 14 






-1-0 39 


— 0^0 


9057 Q3 


Fp TTT 


—0 40 


—0 93 


—0 5^ 


00 




—0 90'^ 


—0 03 


—0 79 


9058 90 


Fp TTT 


—0 38 


— /'i 


—0 5^ 


— 97 

(7. t-> / 




—0 90^1 


00 


—0 QO 


2058.55 


Fe Ill+Cr III 


+0.10 


0.00 


-0.1 


+0.50 




+0.70 


+0.40 


-0.65 


2059.67 


Fe III+Mn III 


0.00 


0.00 


-0.1 


+0.10 


0.0 


+0.70 


+0.45 


-0.99 


90fi3 40 


Mr\ TTT 




—0 18 




—0 01 





—0 90 


-1-0 05 




2065.26 


Fe IIIx4+Fe V 




-0.30 




-0.10 




-0.40 


-0.10 


-0.85 


2066.40 


Mn IIIx2+Ni III 


-0.30 


-0.40 


-0.2 


-0.20 




-0.35 


+0.05 


-0.90 


2068.25 


Fe III 


+0.60 


+0.10 


+0.2 


+0.30 


+0.5 


+0.43 


+0.47 


-1.00 


2069.82 


Fe III+Mn III 


+0.05 


-0.20 


0.00 


-0.4 




-0.15 


+0.20 




2070.56 


Fe III 


+0.47 


0.00 


-0.3 


+0.44 


-0.1' 


+0.89 


+0.55 


-0.79 


2070.98 


Fe III 


-0.35 


-0.35 


-0.5 


-0.19 


-0.1^ 


-0.30 


-0.12 


-0.70 


2074.23 


Fe III 


-0.15 


-0.70 




+0.12 




-0.35 


+0.14 




2076.31 


Fe III 




-0.70 




-0.09 




-0.40 


+0.05 




AVG 




-0.16 


-0.15 


-0.2 


0.02 


0.1 


-0.16 


0.16 


-0.88 


1 a 




0.17 


0.15 


0.2 


0.12 


0.1 


0.23 


0.14 


0.10 



Note. — Abundances are determined from spectrum syntheses using the atmospheric parameters in 
Tables ^ and ^. Dominant features only are identified here. Results that are >2cr from the mean are 
italicized and not included in the average. We find a mean iron-group abundance of [Fe/Hj=— 0.06 ±0.15 for 
the main sequence B-type stars in this Table. This is expected for solar neighborhood objects, and similar to 
the B-type star results in Cunha & Lambert (1994) and Gies & Lambert (1992). 

a,b,c,d,e,f,g,h'p]^g abundances noted are from a blend of more than one feature listed in the table, e.g., all 
entries marked by the letter 'a' are from a single fit to a blend of those features. 

'The Ekberg (1993) wavelength does not fit the feature in PC 0832+676, although it does fit the features 
in the other sharp lined stars. 
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Table 7. Representative Iron-group Abundance Uncertainties 





BD+56 576 


HD 34078 


HD 36591 


HD 216916 


PG0832 




A[M/H] 


A[M/H] 


A[M/H] 


A[M/H] 


A[M/H] 


Areff=±750 K 


TO.01 


±0.2 


±0.10 


±0.01 


TO.04 


Alog g=±0.1 


±0.02 


^0.1 


^0.01 


±0.03 


±0.03 


A,f = ±1 km 


TO.27 


^0.1 


=F0.23 


^0.18 


=F0.23 


A^Ma = ±2 km s-i 


±0.09 


±0.0 


±0.11 


±0.11 


±0.14 


Shift Continuum ±1% 


TO.07 


TO.1 


TO.07 


TO.08 


TO.07 
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Table 8. Boron Abundance Uncertainties 





BD+56 576 


HD 34078 


HD 36591 


HD 216916 


PG0832 




Alog(S/H) 


Alog(B/H) 


Alog(B/H) 


Alog(B/H) 


Alog(S/H) 


Areff=±750 K 


±0.01 


±0.2 


±0.05 


±0.00 


to.oi 


Alog g=±0.1 


±0.03 


TO.1 


±0.03 


±0.04 


±0.05 


= ±1 km 


TO.03 


±0.0 


±0.00 


±0.02 


±0.00 


A^Afa = ±2 km 


±0.02 


±0.0 


±0.00 


±0.04 


±0.01 


Shift Continuum ±1% 


=F0.03 


^0.1 


:F0.06 


^0.03 


=F0.13 


iiB =2 10b 


-0.02 


0.0 


-0.06 


-0.12 


-0.10 


AVrad = ±2 km s-^ 


±0.10 


±0.0 


±0.07 


±0.03 


±0.04 


(AVrad Mn IIIA2065.9)* 


{TO. 10) 


(±0.0) 


{TO. 13) 


{tO.26) 


{to.oo) 



Note. — AB was determined by scaling the Mn III A2065.9 line abundance a priori according 
to the results listed for the five uncertainties examined for the iron-group elements in Table 0. We 
also examine a smaller iiB/^^B ratio since this may be affected by rotational mixing (discussed in 
Section 4). 

'Changing the radial velocity required a completely new fit, thus the new "best" Mn 111 2065.9 A 
abundance is also listed here {in italics). 
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Table 9. NLTE Boron and CNO Abundances in B-type Stars 



Star 





Liter 


ature 





REFS 





Corr 


ected 









B 


C 


N 









B 


C 


N 







THIS PAPER: 
























BD+56 576 




7.89 


7.35 


8.57 




V+00 


2.25^ 


7.84 


7.62 


8.34 


22500 


HD 16582 


<1.3 


8.14 


8.12 




PQOl, 


GL9z 


1.16^ 


8.15 


8.10 




22942 


HD 36591 




8.24 


7.64 


8.64 




GL9/ 


<1.32^ 


8.25 


7.69 


8.67 


26449 


— 


<2.0 


8.25 


7.58 


8.60 


PQOl, 


CL94 




8.32 


7.64 


8.54 


26330 


HD 50707 


<1.7 


8.20 


8.23 


8.87 


PQOl, 


GL92 


<1.5'' 


8.21 


8.29 


8.89 


26768 


HD 205021 


<1.2 


7.96 


7.97 


8.77 




r^T no 


<0.9^ 


7.98 


8.00 


8.81 


25831 


HD 216916 


2.4 


8.15 


7.66 


8.64 


PQOl, 


(jL9z 


2.3P 


8.17 


7.64 


8.61 


23232 


GHRS BORON: 
























HD 886 


2.2 


8.41 


7.78 




PQOl, 


GL92 


2.2 


8.40 


7.75 




21899 




2.29'^ 










P+99 


2.28*^ 








22670 


HD3360 


<1.2 


8.43 


8.26 




PQOl, 


GL92 


<1.2 


8.40 


8.22 




21426 




0.92*^ 










P+99 


0.91^ 








22180 


HD 35299 


2.9 


8.17 


7.71 


8.63 


PQOl, 


GL92 


2.9 


8.19 


7.70 


8.64 


23831 




2.7 


8.38 


7.74 


8.57 


L+00, 


CL94 


2.7 


8.40 


7.67 


8.37 


24000 




2.46'^ 










P+99 


2.41^ 








24670 




2.84*= 










C+97 


2.84^^ 








24000 


HD 35039 


2.6 


8.40 


7.76 




PQOl, 


GL92 


2.6 


8.35 


7.70 




20547 




2.92^= 


8.50 


7.85 


8.60 


C+97, 


CL94 


2.92^ 


8.36 


7.65 


8.34 


20550 


HD 36285 


2.43^= 


8.57 


7.95 


8.80 


C+97, 


CL94 


2.43^ 


8.48 


7.77 


8.55 


21930 




1.8 










PQOl 


1.8 








2223ff^ 


HD 36430 


2.53^= 


8.54 


7.89 


8.84 


C+97, 


CL94 


2.53^ 


8.38 


7.67 


8.57 


19640 




2.5 










PQOl 


2.5 








19640 
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Table 9 — Continued 



Star 






Liter 


ature 





REFS 





Corr 


ected 









B 




C 


N 









B 


C 


N 
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Note. — This table summarizes the NLTE B, C, N, and O data for B-type main 
sequence stars, excluding stars with 18,000 < T^s < 29,000 K (see text). The left 
side of the table quotes NLTE abundances as published in the literature. The right 
side of the table corrects the NLTE CNO abundances by (i) including adjustments 
to correct for the use of Gold, instead of Kurucz, model atmospheres, and (ii) 
reducing the GL92 temperatures by 3.4%, and adjusting their abundances according 
to the offsets listed in their Table 9. This reduction in the GL92 temperatures has 
a negligible effect on the B III abundances (e.g., see Table There are some 
significant differences between the GL92 and CL94 abundances for the five stars in 
common. The differences in oxygen appear to be due to differences in ^, whereas 
those in C appear when two strong C III lines are included by GL92. N is in good 
agreement throughout, with the only significant difference of ~0.1 dex occuring for 
the hot star HD 36960. 

'^Boron abundances from this paper. 

''Proffitt et al. (1999) LTE B III abundances have been NLTE corrected 
(corrections are from this paper, see Table ^) for this table. On the right side 
of the table, their abundances are also corrected for the differences between their 
temperatures and the reduced GL92 temperatures listed (small effect). 

'^Boron from B II 1362.5 A, which is severely blended with metal lines. 

'^PQOl temperatures are the same as CL94 or they have adopted GL92— 820 K. 
This offset is very similar to our reduction of 3.4%. We note when the PQOl and 
GL92 temperatures differ by more than 250 K. 
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Fig. 1. — Coadded i^STSTIS spectra for three B-type stars (thick line) and their spectrum syntheses 
(thin line). The observed spectra were smoothed for a 3 pixel resolution element, and the higher 
resolution from the E230H grating can be seen for HD 16582. The iron-group metallicities in Table ^ 
were used for each synthesis. Features listed in Table ^ are identified. Interstellar lines are marked; 
PG 0832+676 has two sets of IS lines. 



ts5 



C/3 




2060 

Wavelength (A) 



p3 



C/3 




00 
C5 



2070 2075 
Wavelength (A) 



o 



O ^ 03 

00 tr- 
00 S P 



+ ^ 

fD 



95 



o 

95 



I— I 03 
t/1 



CD 
CC 
CD 
U2 

CD 

i-i 

CD 
CD 



CD 



ora CD 



td 



H- 
o 



CD 



o 
o 



95 

go' 
O 



CD 
r+ 
CD 

n 

CD 

a 

95 

o 

o 

00 

w 
to 

+ 

-a 



1.0 



1.0 



w 
c 

CD 



.> 1.0 

CD 

DC 



1.0 



0.0 



"1 r 



"1 r 



"1 r 



"1 r 



h'd 216916 




PG 0832 




J L 



J L 



J L 



J L 



2065 



2065.5 



2066 2066.5 

Wavelength (A) 



2067 



2067.5 



2068 



Q M ^ 



o 
n 




00 

Or • 

CO 



cr 

O CD 
cc' 

JO ^ 

a- 

o 



CD 



!=r o 

o B 



td 

o 

o 

1:3 

'< 



3^ 



3? S 



o tr 



O 



CD 



CD 



CO 



CD 



CD 
O 



05 

fl- 



ee 
CD 

cu 



CD 



0} 
> 

DC 



1.0 



0.0 



"1 r 



"1 r 



"1 r 



"1 r 



HD 34078 




. . . 



HD 205021 





HD 50707 




4 




J L 



J L 



J L 



J L 



2065 



2065.5 



2066 2066.5 

Wavelength (A) 



2067 



2067.5 



2068 




Fig. 6. — Boron syntheses for HD 36591; 12+log(B/H)LTE = l-05, 1.45, 1.85 are shown. Boron is 
clearly weak in HD 36591. For comparison, an HST GHRS spectrum for HD 35299 is shown. These 
two stars are both in the Orion OBI association and have very similar atmospheric parameters, 
yet Proffitt et al. (1999) and Lemke et al. (2000) fit a normal LTE boron abundance to HD 35299 
(12+log(B/H)=2.5 and 2.7, respectively). Interstellar lines are marked for HD 36591. 
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Fig. 7. — Predicted equivalent widths for boron lines in B-type stars, from NLTE calculations. The 
B III 2065.8 line strength plateaus between 18,000 K and 29,000 K. At hotter temperatures, an 
increase in T^g of only 5 % reduces the predicted equivalent width by 50 %, or 0.4 dex in the boron 
abundance. 
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Fig. 8. — Boron abundances versus temperature. Temperatures from Gies &: Lambert (1992) 
have been reduced by 3.4% (see text). This plot shows that there is no significant trend in the 
B III 2065.8 A abundances, although we note the upper and lower temperature limits where the 
B III 2065.8 A line yields the most reliable abundances. Thick lines connect boron abundances 
found for the same stars from different analyses. 
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Fig. 9. — NLTE corrections to UV boron lines. These calculations arc identical to those discussed 
in Cunha et al. (1997); here, we extended the grid to hotter temperatures (35,000 K) for the B III 
features. The change in orientation of the B III NLTE correction in the hotter models appears to 
be due to the much weakened boron lines being dominated by overionization effects (as opposed to 
line transition rates, which dominate at lower temperatures). 
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Fig. 11. — Gravity versus T^s for B-type main sequence stars in Table Open diamonds mark 
normal stars (B>2.2 and N<7.8), filled diamonds show boron depleted stars (B<2.2 and N<7.8), 
and asterisks denote the N-rich stars (B<2.2 and N>7.8). Evolution tracks are from Heger & Langer 
(2000) for stars of 8, 10, 12, 15, and 20 Mq, two rotation rates, and considering ^-barrier effects, 
through the H-core burning phase. This plot shows the mass and age range of the B-type stars 
examined in this paper, as well as the distribution of boron depletions and nitrogen enrichments 
on the main sequence. 
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Fig. 12. — Predictions of boron versus nitrogen (left panel) and boron versus main sequence age 
(right panel) from the rotating stellar evolution models from Heger & Langer (2000). Predictions 
for a 12 Mq model with five rotational velocities are shown {dotted lines, velocities are marked in 
the right panel). For comparison, predictions for a 15 Mq model with rotational velocities of 200 
and 450 km s~-^ are shown {solid lines). Additionally, 10 and 15 Mq models, rotating at 200 km s~^, 
and that ignore the effects of ^u-barriers arc shown. In the left panel, it is striking that the boron 
versus nitrogen relationship is nearly the same for various masses, rotation rates, and efficiency of 
mixing with/without ^-barriers, although some tracks end before significant abundance changes. 
However, in the right panel, it can be seen that the rate of boron depletion is sensitive to these 
parameters. 
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Fig. 13. — Boron versus nitrogen in main sequence B-type stars, with predictions for the 12 Mq 
models at 200 km s~^ through the H-core burning phase, and three sets of initial abundances [(Bj, 
Ni) = (2.6, 7.6), (2.8, 7.7), and (2.8, 7.8), all within range of the interstellar and ZAMS stellar 
abundances]. The stellar abundances plotted are those listed in Table |^. Thick lines connect boron 
abundances for the same stars from different analyses. The observations are in good agreement 
with the predictions, suggesting that rotational mixing is significant on the main sequence in some 
B-type stars. 
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Fig. 14. — Boron versus cluster ages (see Table ^ for ages and references) for B-type stars in 
associations. Stellar data symbols are the same as in Fig. 13. Tracks are the same as in Fig. ^ (some 
tracks are not shown here for clarity). One additional track is included to show the predictions for 
a 20 Mq star with ^u-barrier effects, rotating at 450 km {solid line). Three boron depleted young 
stars (HD 36591, HD 50707, HD 205021) are fit well by the rotating model predictions, although 
the most boron depleted stars are best fit by the rapidly-rotating 20 Mq model, which is likely too 



massive for these stars (see Fig. |ll| ). Many of the older stars in this sample that do not show boron 
depletions are less massive (e.g., the 10 Mq track, dashed line, fits them well), though some are 
likely to be true slow rotators as well. 



